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ABSTRACT 


An investigation of the intensity fluctuations of 28 pulsars 
near 0.4 GHz indicates that scintillation spectra have a Gaussian 
shape, scintillation indices are near unity and the scintillation 
bandwidth depends linearly on dispersion measure. Observations 
near 2.5 GHz suggest a strong dependence of the frequency at which 
scintillation indices fall below unity on dispersion measure. Recent 
multistation measurements of scintillation provide values 
or limits for the scale size of the scattering diffraction pattern. 
The dependences of scattering parameters on dispersion measure is 
discussed in terms of the current models.. It is suggested that any 
line of sight through the galaxy encounters increasingly -rare , 
increasingly- large deviations of thermal electron density on the 
scale of 10 11 cm. 
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I . INTRODUCTION 


Scheuer (1968) and Salpeter (1969) interpreted the rapid, 
narrow-bandwidth , random fluctuations of pulsar signals first 
observed by Hewish and his colleagues (19 68) a.s the result of 
scattering by minute density irregularities in the column of 
interstellar thermal electrons that was known at the time to 
disperse the pulsed signals. Rickett (1970) extended the model 
o f Scheuer using the work by Ucsinski (19 69) and mad e a detail ed 
comparison of his observations to the theory. Mutel et al„ (1974) 
have recently discussed power- law density irregularity spectra 
in the interstellar medium based on the work of Salpeter and 
Lovelace (1970). Cohen and Cronyn (1974) describe observable 
consequences of interstellar scattering based on current data and 
theory. In the models the thermal plasma is assumed to be uniformly 
distributed on a galactic scale, with rare exceptions, since HI 
absorption measurements for pulsars provide distance estimates which, 
with the electron column densities (dispersion measures DM), yield com- 
parable average electron densities along different lines of sight (Gordon 
and Gordon 1973, Gomez-Gonzalez and Gueiin 1974). Also, pulsars lie below 
a line DM «= csc(jb|) in the -DM galactic latitude(b) plane (Hewish 1970, 
Condon and Backer 1974), and above a line L « DM in the 
Luminosity(L) -DM plane (Cornelia 1973); both of these results lend 
support to the idea that the dispersing plasma is distributed generally 
in the interstellar medium. The distribution of irregularities within 
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the dispersing plasma is also probably general in 

the sense of frequent, ’on the scale of 10 pc, perturbations of 

the electron density as discussed by Salpeter, rather than in the 

sense of isotropic turbulence. Scheuer and Salpeter considered 

thin-screen approximations to interstellar scattering, while 

Rickett applied the work of Uc inski which deals with extended screen 

effects . 

The rare exceptions mentioned above are the thermal plasmas 
in the Crab nebula, in the Gum nebula, and in miscellaneous HI I 
regions (Prentice and ter Rear 1969). Although Gomez-Gonzalez and 
Guelin (1974) do not find a dependence of 'the mean electron density 
on distance, the average perturbation of the density about the mean 
could increase with distance as the line of sight intersects 
increasingly-rare , increasingly- large deviations of the electron 
density, particularly in directions toward the galactic center. 

The frequency dependence of scattering parameters (Rickett 
1970, Lang 19 71, Rankin et _aj_, 1970, Komesaroff et ai. 1970, Backer 
1974, Mutcl et ajl. 1974) is consistent both with Scheuer's model 
which assumes a homogeneous and isoLropic medium and a Gaussian 
irregularity spectrum and with Salpeter 1 s model which includes the 
effects of a density irregularity spectrum with a logarithmic s^ope 
of S = 4 (see Mutcl et al.)» The frequency dependence of the 
ex tended -screen model (Rickett 19/0) is also in agreement with the 
observations. The dependence of the scattering parameters on 
distance, or dispersion measure for a uniform interstellar medium, 
is less clear. Only one parameter, the correlation bandwidth (B s ) > 
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or its Fourier transform conjugate the pulse broadening time, is 
clearly dependent on dispersion measure (Rickett 1970, Sutton 1971). 
Rickett finds a relation of B * DM -2 which is consistent with the 


s 

above models, while Sutton (1971) finds a trend of “ DM - ^. In 
this paper we investigate the dependence of other scattering 
parameters on dispersion measure: the correlation time (scintilla- 

tion bandwidth), the diffraction pattern spatial scale, and the 
critical frequency for scintillation. Before presenting these we 
discuss the shape of scintillation spectra. Additional topics are 
given: the use of scintillation of malnpulse and interpulsc 

radiation as a means of resolving the transverse separation of the 
two emitting regions, and a simulation of scintillation using a computer 
code developed by R. Lovelace (1974). Finally we discuss the observed 
strong dependence of scattering parameters on dispersion measure 
and present a model scintillation spectrum with explicit dependence 
on DM and radio frequency. 


II. OBSERVATIONS 

Pulsar recordings extending over intervals ranging from an hour 
to several days were made over the past few years at several 
observatories (National Radio Astronomy Observatory 1 , Jodrell Bank 
and Arecibo Observatory) to investigate pulse-to*pulse fluctuations, 
flux-density spectra, interstellar scattering diffraction patterns, 
individual pulse polarization phenomena and differential positions. 
Table 1 summarizes these experiments. This paper presents the salient 


^NRAO is operated by Associated Universities, Inc. , under contract 
with the U.S, National Science Foundation. 
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bostures of interstellar scintillation contained in this large 
collection of data, Other papers have discussed or will discuss 
the individual experiments with more detail, although relevant 
aspects of some observations not given in table 1 will be presented 
in the succeeding sections of this paper. All the data described 
in table 1 provided the author with experience on the fluctuations 
of puxsar signals, only part of the data will be discussed here due 
eithei to the redundancy of the information, to the incompatability 
of magnetic tape formats, or to the fact that the radiometer bandwidths 
or the lengths of recording were insufficient for the analysis of 
scintillation, 

III * SHAPE OF SCINTILLATION SPECTRUM 

The shape of the scintillation spectrum, which gives the variance 
density of scintillation as a function of fluctuation frequency for a 
pulsar, is confused at some level by the variability of the emitted radia- 
tion, due to changes of parameters intrinsic either to the pulse emission 
mechanism or to propagation conditions in the near pulsar environmei 
Variance-density spectra of pulsar data typically have the following 
constituents (Lovelace and Craft 1968, Backer 1973) which are labeled 
in a schematic spectrum in figure 1: 

(a) a rapid decrease near zero, followed by 

(b) a slower decrease out to a point which varies from 0.1 
to 0,5 cycles per pulse period (c/P^), and beyond which 

(c) a constant, or "white", level is found with 
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(d) one or more quas iperiod ic features superposed on the 

above "continuum", ‘ 

Additional increases in pulsar spectra below frequencies of inverse 
weeks (2 x 10" Hz), which have been observed at meter wavelengths 
(Cole, Hesse and Page 1970, Huguenin _et ajL 1973), will not be 
discussed here* Constituents (c) and (d) are properties c£ the 
pulsar since their corresponding temporal fluctuations are correlated 
over broad bandwidths, are related to complex changes in total-intensity 
and polarization emission waveforms unlike the random variations 
expected from scintillation, and .are in case (d) periodic, Rickett 
(1970) and Lang (1971) showed that constituent (a) is dependent on 
dispersion measure(DM) , and therefore to be- identified with scintilla- 
tion arising in the interstellar medium: they found that the fluctua- 

tions producing (a) are reduced considerably by observations with a 
radiometer bandwidth exceeding a characteristic correlation bandwidth 

B which varies with DM^ ^ (see also Sutton, 1971) < 
s 

Since finite lengths of data are used to estimate pulsar fluctuation 
spectra, the sampling function could lead to a spurious component { b) 
due to the excess modulation at low frequencies. This is judged not 
to be a serious problem based both on observations of constituent (b) 
when scintillation has been removed by the use of large radiometer 
bandwidths and on computer simulations of scintillation records which 
are discussed in tVIII.. • Often constituent (b) can be identified with 
intrinsic variability when the fluctuations are related to deep and 
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abrupt nulls (Backer 1970a), equally abrupt mode changes (Backer, 

1970b, Lyne 1971), isolated bursts of emission (Huguenin, Taylor 
and Troland 1970, Backer, Rankin and Campbell 1974),. and low- 
frequency periodicities (bang 1969), most of which are correlated 
over hundreds of MHz. This constituent usually follows an exponential 
decrease from zero frequency, and is often sufficiently strong to 
confuse measurements of scintillation time scales. Thus a careful 
inspection of the low-frequency end of a spectrum is needed to 
determine the shape of the scintillation spectrum. One can be 
guided by the hypotheses that (i) scintillation spectra will be 
Gaussian for interstellar electron-density spectra which are consistent 
with a radio-frequency ( v) dependence of the correlation, bandwidth B > 
or the inverse of the pulse-broadening time, of (Rickett 1970, Rankin 
et al « 19 70, Komesaroff ejt al_. 1970, Lang 1971, Backer 1974), and of the 
scattering angle of (Mutel et al. 1974, Vandenberg 1974); ( ii) 
normalized spectra should have unit area corresponding to strong 
scattering in the Fraunhofer, or far field, limit; and (iii) any 
consistent measure of the width of the spectra should decrease with 
v and increase with DK V , where a is 0,5 for constant motion of the 
observer relative, to a uniform medium of Gaussian irregularities 
(Scheue.r 1968) or for the power-law irregularity spectra discussed 
by Salpeter (1969) with S = 4 as required by the above dependences 
on radio frequency, 

A study of the spectra of 28 pulsars at radio frequencies near 
400 MHz, which is presented in 'j IV, reveals that constituent (a) 
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exists in all but three cases for which either the length of the 
recordings was too short or the radiometer bandwidth was too large 
to observe scintillations. Furthermore, in many cases the exponential 
constituents are weak enough to fit Gaussians to the scintillation 
spectra over several octaves of fluctuation frequency, from a few 
percent of the e width, f £ , out to t e * 

In figures 2 to 4 variance-density spectra are given for three 
pulsars. In these plots both high and low resolution spectra are 
given with Gaussian curves (dash-dotted) , selected by eye from a 
template of curves, drawn through the low frequency ends of the high 
resolution spectra (solid). The low-resolution spectra (dashed) provide 
comparison levels, the levels of uniform, or "white", fluctuations 
at, or near, the Nyquist limit* Error bars indicate the uncertainties 
expected due to limited smoothing of adjacent points in the spectra... 

For PSR 0823+26 (f ig. 2) the high-resolution spectrum from 408 -MHz 
data decreases to the instrumental noise level, rather than to the 
level of the white fluctuations given by the companion low-resolution 
spectrum based on 4 30 -MHz data. In the other two cases instrumental 
noise is insignificant. These variance-density spectra are defined 
in an absolute manner to allow comparison at other radio frequencies 
or other epochs. Comparison functions are given in figure 5 from 
2 69 5 -MHz data recorded in 1973 May (see table 1). For PSR 0329+54 
the width of the scintillation portion of the spectrum is scaled 
down from the 408 -MHz result in proportion to the ratio of the radio 
frequencies as expected; the area of the scintillation spectrum 


8 



corresponds to a scintillation index of only 0.4 which is small 
compared to the expected value of unity. Scintillation at 2695 MHz 
is discussed in more detail in 2 V , and scintillation index variability 
is discussed in 9 VIII. The level beyond the scintillation portion, 
f > 0.3 mHz, is higher than that expected based on the 408-MHz 
data above 3 mHz (fig. 2) which suggests either an exponential con- 
stituent of fluctuations or a radio-frequency dependence of the continuum 
level. The 2 69 5 -MHz variance-density spectrum for PSR 0823+26 has a 
scintillation portion which agrees with the 408-MHz data (fig. 3) based 
on the conventional notions of frequency scaling described above. The 
level beyond the scintillation portion is undetermined due to instrumental 
noise and the low Nyquist frequency. The 2695-MHz spectrum for 
PSR 1133+16 has a scintillation portion which also agrees with the 408-MHz 
data, if account is taken of both weak scattering (see $V) and 
the level and slope of the portion of the exponential constituent 
of the 408-MHz spectrum just beyond 3.5 mHz (fig. 3). 

Further support- to the argument that the exponential constituents 
of pulsar variance-density spectra are due to "intrinsic” causes is 
that neither the transition point where the Gaussian scintillation 
spectra and the exponential constituent are equal, nor the absolute 
level of the exponential constituent hear any relationship to 
dispersion measure Tor the sample of 28 objects discussed here. In 
contrast, the width of the Gaussian portion does depend on dispersion 
measure as shown in the following section* 
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It is concluded that all pulsar data are consistent with a 
Gaussian shape and unit area for the normalized, variance-density 
spectra of interstellar scintillation. There is no evidence for 
the nongaussian behavior discussed by Little and Matheson (1973) 
based on the work of Lang (1971). 

IV. WIDTH OF SCINTILLATION SPECTRA 

Rickett (1970) measured scintillation time scales, t , for 
12 pulsars from the separation of peaks on analog recordings and 
found, roughly, an inverse relation between t and DM, There are 
several reasons not to expect a rigid relationship: 

(a) insufficient data for a stationary result; 

(b) nonuniform medium of irregularities in electron density; and 

(c) variation from object to object of the transverse velocity 
(or velocities?) of the scattering medium with respect to 
the line of sight to the pulsar. 

The recordings near 400 MHz listed in table 1 were analyzed to 
identify the scintillation spectra, as discussed in the preceeding 
section, and then to measure the areas, and the widths at 

£ , of the Gaussian scintillation spectra,, The results are summarized 
in table 2 which gives systematic pulsar names (PSR) , periods (P) , 
dispersion measures (DM), lengths of recording in periods (N x M, where 
N is the number of samples and M is the sampling interval), square 
roots of the areas of the scintillation spectra (m, the scintillation 
index), e widths of the scintillation spectra ( f ) , transition 
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frequencies where the scintillation spectrum and the intrinsic 
spectrum are equal (£ t ) , sources of the data by observatory (AO = 
Arecibo Observatory, GB = Green Bank, JB - Jodrell Bank) and by 
experiment (IF — fluctuation , IS ~ scattering , FN — polarization) , 
and the qualities of the measurements (f x N x M x P = f e x T) . 

Gopies of the spectra may be obtained from the author upon 
request. 

The values of from table 2 are plotted against dispersion 
measure in figure 6 0 Only upper i units were obtained for 
PSR's 1541+09, 1604-00 and 1944+17 due either to the brevity of the 
recordings or to the narrowness of the correlation bandwidth relative, 
to the radiometer bandwidth. A value of f e for PSR 0833-45 
(DM ~ 70 pc era' 3 ) at 40 8 -MHz is included in figure 6 based on 
work at 834 MHz (Backer 1974). 

The trend of the 400 -MHz data is given by 

( 1 ) = 0 S 04 DM v _i MHz, 

where DM is given in units of pc cm ^ and v in GHz here and in 'follow- 
ing equations; individual points are spread by a factor of three 
from this trend due either to large scale (~ 100 pc) inhomogeneities 
of the scattering material or to discordant velocities as discussed 
at the beginning of this section. This strong dependence of f g on 
DM evident in figure 6 is not expected within the context of uniform- 
medium models of scattering. Further discussion may be found in SIX. 
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V. THE CRITICAL FREQUENCY FOR SCINTILLATIONS 


Observations oi scintillation at mater wavelengths are consistent 
with scintillation indices of unity (see i^IV and references therein). 
Indices of one are expected when ( i) random-phase deviations imposed 
on a plane wave by the interstellar electron-density irregularities are 
many radians, and (ii) the observer is sufficiently far from the 
density, or phase, irregularities to allow the reception of scattered 
radiation from many paths. At some radio frequency for any object, 
the critical frequency, one of these conditions will be violated, 
and the scintillation index will drop below unity. For a collection 
of objects, at a given radio frequency, the index will fall below 
unity for objects closer than some critical distance, or dispersion 
measure, assuming a uniform medium. Explicit statements of the above 
conditions involve the product DlT with <Xj = 1/2 and = -1 f 0r 
condition ( i) and = 1/2 and $2 - "2 for condition (ii). Thus, 
given the critical dispersion measure for weak scintillation DM_ at 
some radio frequency v, one can obtain the critical frequency v as 
a function of DM for the above models. 

Downs and Reichley (1971) were first to show that at least one 
of the conditions above was violated for some objects at 2385 MH z. 

Their scintillation indices are upper limits for interstellar 
scintillation, since they did not consider contributions from intrinsic 
effects ( § III) ^ Their data are given as published in figure 7(f). 
Indices have also been estimated from the scint il iat ion portions of 
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variance-density spectra of the 2 69 5 -MHz data listed in table 1. 

The results are listed in table 3 and plotted in figure 7 for the 
1973 May flux-density experiment (A arid ® ) and the 1974 May position 
experiment (B and A) „ The latter data sets extended over longer periods 
of time than the former and therefore are probably more reliable when 
available* For PSR’s 0833-43, 1749-28, and 19334-16 the correlation 
bandwidths are known (Sutton 1971, Backer 1974) and corrections to the 
indices for bandwidth suppression, have been applied (O'), 

Below the critical dispersion measure, at a given radio frequency 
such as 2695 MHz, the scintillation index will decrease in proport ion 
to DM 11 where a is 1/2 H condition ( i) is violated and 1 if both 

conditions are violated, or if condition (ii) is violated sufficiently 
strongly (Lovelace 1970). A simple model of this cutoff is 

1 

(2) - “ a. 

1+ (DM c /DM) 

The data is not sufficient to determine within the range of 0.5 to 1.0 
discussed above , although it appears that a value closer to 1 than to 
0*5 is favored with DM = 5 pc cm (sc.e fig. 7). Taking n. = 1 and 
p = 1 (arbitrarily) , one finds that the critical radio frequency for 
scintillations will vary as 

(3) v c = 0.5 DM GHz, 

although both the constant and the dependence on DM are clearly uncertain 
from the above discussion. 
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In 1974 May 53 objects were observed at 775, 968 and 1440 MHz over 
an eight “"day interval (table 1). At each frequency three or four 
measurements of pulse energy, averaged over about four minutes, were 
obtained on as many separate days. For objects with dispersion 
measures less than —25 pc cm 3 the correlation bandwidths of 
scintillation exceed the radiometer bandwidth of 20 MHz and therefore 
the samples are independent realizations of scintillation plu s day- 
to-day intrinsic variations. For objects with dispersion measures 
exceeding —80 pc cnf 3 , the scintillations are quenched by the radio- 
meter bandwidth and the samples exhibit only the day-to-day varia- 
tions. (The data for objects with DM -80 pc cm -3 confirm the steep 
dependence of correlation bandwidth on DM discussed by Sutton (1971)). 
Since in the latter case (DM > 80 pc cm 3 ) the rms deviations about the 
raaan pulse energy are less than 0.2 times the mean, we conclude that 
the rms deviations of the former group, DM < 25 pc cm ' , are dominated 
by interstellar scintillation. Figure 8 gives the radio frequency 
dependences of the rms deviations for nine objects with DM < 13 pc 
from the 775-1440 MHz experiment and from the data in tables 2 and 3 
at 400 MHz and 2695 MHz, respectively. This figure confirms the 
existence of a critical frequency for scintillations and suggests 
that the constant in equation (3) should be closer to 0.1. 

Equation (3) is probably accurate enough to indicate that it will be 
difficult to study many pulsars in the weak scintillation regime. 
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VI. SCATTERING DIFFRACTION PATTERN SCALE 

The intensity diffraction pattern produced in the plane of the 
observer by interstellar scattering below the critical frequency of a 
point source of unit strength may be described by its wave number spectrum 


(4) 


M(q) = exp[-a 2 (q x 2 + q 2 )/2], 


where a is the scale size of the diffraction pattern and q = (q„,q f ) 
is the vector wave number. The gauss ian shape of M(q) is inferred 
from the gauss ian shapes of the variance "dens ity spectra D( £) described in 
III under the simple assumption that a scintillation time and a pattern 
wav i- numb er are related by a speed. With this ass ump t io n , t he follow- 
ing relations hold: 


D(£) D( ») 

A -x 


(5) 


M( q) 


M( r) , 


where refers to Fourier tranformat ion and » refers to a one- 

dimensional integration along the coordinate perpendicular to the pattern 

velocity u and to a change of variables 2nf/ju| = q. The size a varies 
- 1/2 

as DM v for the uniform-medium models of Scheuer (1968) and Salpeter 
( 19 69 ) «. Multiple station measurements of puls ar sc int i 1 1 at ion are 
consistent with values of M(rJ which are less than unity: Slee et al . 

(19 74) observed seven objects at 32 7 MHz with |,r| - 8000 km, and the 
unpublished results of Backer, Lyne and Zeissig at 408 MHz with 
ixl x' 5300 km (table i) are consistent with those of Slee et al. We 
estimate a front these observations to be 
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( 6 ) 


a = (1.4 ^'5) ^ loi ° DM" 1/2 cm, 
although, given the results of §§IV and V, a stronger dependence on DM 
might be expected. Equation (6) is valid for v < > when v '> v c one may 

substitute v = v into equation (6). The uncertainty in a is based on 
the paucity of measurements and, presumably, on inhomogeneities of the 
scattering. 

VII. SCINTILLATION OF MAINPULSE AND INTERPULSE EMISSION 

Investigations of ionospheric and interplanetary scintillations 
have often been a battle to separate medium effects from soprce structure 
effects since the "resolving power" of these scint i 1 1 at ions is in the same 
regime as the angular structure of many extraterrestrial, natural 
sources. The situation is further confused by the time-variable 
properties of these scattering media. The "resolving power" of. the 
interstellar medium is generally thought to fall in the range where 
incoherent radio sources at the present level of detection are 
relatively large, and therefore do not scintillate, and all pulsar 
emitting regions (and perhaps some molecular -line masers) are 
relatively small, and therefore scintillate 100 per cent below the 
critical frequency ( 9 V). We explore the first point elsewhere (Lovelace 
and Backer 1972 , Condon and Backer 1974 ) and the second point below. 

The resolving power of the interstellar medium compared to an 
angular size scale for a pulsar emission region may be expressed as 
the ratio of twice the scattering diffraction pattern scale ("VI) 
to the light cylinder radius of the pulsar, r c = c P/ 2 n, 

(7) 2 a/r c - 6 DM L ^ v P~ x . 
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This ratio is 0.6 for v = 0.4 GHz, DM = 16 pc cm" 3 , and P - 1 seCo 
Thus sources of pulsar radiation from any one pulsar which are distributed 
over a distance r c at the time of emission to the observer can be 
distinguished by their independent scintillations. One is essentially 
collecting radiation scattered by many irregularities in the intervening 
medium with propagation paths which are phase stable tor one period 
so that differential positions can be obtained for radiation arriving 
at different times in the pulse period in a manner analagous to long- 
baseline interferometry. Most pulsar radiation is confined to a few 
per cent of the period and therefore is probably from a source distributed 
over not more than a few per cent of pulsar latitude and longitude. 

However, several of the short -period objects have interpulse emission 
which is spaced roughly 180 degrees of longitude, as observed, from 
the mainpulse emission region (Backer, Boriakoff and Manchester 1973 
and references therein) „ The two regions then could be displaced in 
1 ati t ude by an angle such as twice the angle be tween the magnetic d ipo le 
axis and the rotational axis, 0. The regions, when emitting to the 
observer, would be displaced by 2 r tan 0, where r is the distance from 
the center of mass of the pulsar to the emitting region, r ^ r c . 

For PSR 0823+26 and PSR 1929+10 there was no evidence for 
independent scintillations of the mainpulse and interpulse emissions 
in the relatively short recordings of the Arecibo polarization survey 
at 430 MHz (table 1). This observation sets limits of 1.2 and 8.0 r c , 
respectively, on the separation of the emitting regions at the time of 
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emission. These limits are not interesting for discussions of pulsar 
emission. Indeed, one can only improve them (i) by observing at lower 
radio frequencies where sensitivity would be a serious problem, or (ii) 
by observing higher dispersion measure objects, such as PSR 1055-51 
(McCulloch et al. 1974) which has a strong interpulse, but a short period 
(eq, (7)), or (iii) by finding longer-period objects with interpulses such 
as the elusive (or dead?) PSR 0904+77 which may have an interpulse and a 
period of 1.5 s, but which has not been detected since the months of its 
discovery (Taylor andHuguenin 1969). 

The problem can be turned around in the case of a short -period binary 
pulsar (Taylor and Hulse 1974) . In this case we know the path of the 
radiation source and can probe in principal the stability of the dispersing 
and scattering medium during an orbit. 

VIII. COMPUTER SIMULATION OF SCINTILLATION 

R. Lovelace (1974) has developed a computer code to simulate 
far-field scattering from a one-dimensional thin screen with an arbitrary 
spectrum of phase fluctuations. A version of this code has been developed 
to look at the intensity fluctuations from a scattering screen with a 
gaussian distribution of phase fluctuations whose root-mean-square (rms) 
deviation from zero is 20 radians and whose scale size is 10 11 cm at a 
distance of 500 pc and at a radio frequency of 0.4 GHz. The intensity 
data in the observer's plane resulting from illuminating this screen 
with a plane wave was broken into sections, in an attempt to mimic 
actual observations with f g T ~ 5 as discussed in §§III and IV. In 
each section a scintillation index and bandwidth was determined. We 
found that the scintillation index is less well-defined than the 
scintillation bandwidth in the ensemble of measurements. For the 
above example, with <f e T> = 5 in the ensemble of 16 measurements. 
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estimates of the index had a mean of .95 and an rms deviation from 
the mean of 0,4, while estimates of 'f varied from 5/1 by only 1/T. 

The gaussian shape of the resultant scintillation spectra obtained 
over three orders of magnitude in variance density, or three octaves 
of fluctuation frequency. This supports the statement that constituent 
(b) discussed in $111 is not the result of a smearing of the 
scintillation spectrum by the transform of the sampling function. 

The stability of f in this experiment is reasonable since with a 
gaussian spectrum all "lobes" of scintillation have similar widths 
while their amplitudes have an exponential distribution yielding 
the greater uncertainty in the index. 

IX. DISCUSSION 

The pulsar observations discussed in this paper suggest that 
scintillation indices for point source observations at meter wave- 
lengths are near unity, and that scintillation spectra are gaussian 
in shape and have a width which scales linearly with dispersion 
measure,, The scatter of individual objects from the spec trum wid th( f ) - 
dispersion measure(DM) trend is a factor of three (fig. 6), and is 
attributed both to scatter of the mean electron density along various 
lines of sight (Gbmez-Gonzalez and Guelin 1974) and to nonuniformities 
of the fractional deviation from the mean electron density in different 
regions of the galaxy. Also, the observations suggest a strong 
dependence of scintillation index on dispersion measure in the regime 
of weak scintillation which has been translated to a strong dependence 
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of the critical radio frequency which separates weak and strong 
scintillation regimes on dispersion measure. Sutton (1971) has found 
a similar strong dependence of correlation bandwidth on dispersion 
measure from a compilation Of published data which has been confirmed 
by the measurements discussed in $V. 

In Scheuer's model (1968) and in Salpeter’s model (1969) 
the scintillation bandwidth and correlation bandwidth have the 
following dependences on velocity of scattering medium trans 
verse to the line of sight (v), thickness of the medium (z), rms 
density perturbation of thermal electrons on a scale of 10 cm 
(A) and radio frequency (v) in the strong scattering regime : 


(8a) 


f (X V 

e 


1/2 

z 


v , and 


(8b) ■ B a z -2 A“ 2 , respectively, 

s 

We observe that f o: DM and B cc DM“^„ If DM qc z, statistically, as 
e s 

argued by Gomez -Gonzalez and Guelin and others, then equation (8) 
leads immediately to A « z and v 2 ‘ 1 ^ 2 (if B g <x DM" 3 , then we could 
have A a: z l/2 and v ~ constant)., On the other hand, if DM « z 
then we could have both v and A be independent of DM., and therefore z 
Since the HI measurements favor DM z , we conclude that the former 
relations of A( z) and v(z) are closer to the truth. That is, the 
typical deviation of electron density in the Interstellar scattering 
centers is a growing quantity and the effective pattern speed is 
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either constant or slowly decreasing with increasing path length 
in the interstellar medium. A decrease of speed with z is difficult 
to understand since differential galactic rotation will certainly 
lead to an increasing speed with distance. This puzzling result 
suggests the need for firmer estimates of the dependence of f^ and Bg 
° n DM and for more thought on the expected behavior of these 
quantitives for realistic models of the galaxy. 

The explicit formulae for f g (DM,v) and v c (DM) presented 
in this paper, and the firm assessment of the shape and normalization 
of scintillation spectra allow the construction of a model scintilla- 
tion spectrum for a point source: 

(9) P(f) = (-— r-) 1/2 U + 1 ex P (-f 2 /2£ 2 ), 

nf e c 

where f (DM.,v) is given by equation (1) and v c (DM) is given by 

equation (3). The computer modeling of scintillation discussed in 

9 VIII suggests that a single recording of length T = 5/f will follow 

equation (9) with a normalization (scintillation, index) error of 

< 0*4 and width error of ^ 0.2 f e , where the normalization of unity and 

width of f are ensemble averages of many independent recordings ■ 
e 

The observations upon which this paper are based confirm the computer 

model: when f T is small, there are large deviations of normalization 

e 

and width; when f T is large, there are small deviations . Equation 
e 

(9) can be used to predict the scintillation of a point source of 
known dispersion measure, and, with the diffraction pattern scale 
size (eq. 6), to predict the scintillation of an extended source 
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(Condon and Backer 1974) * For example, if the distances to galactic 
radio stars or molecular line masers are known, then it should be 
possible to estimate a dispersion measure and consequently to predict 
an expected scintillation time scale. 

The effective pattern speed is the product 2rra f , or 

+200 _i 

50 ( ) krns for DM = 16 where the bulk of the measurements exist. 

The above discussion led to a speed dependent on DIP with 0 > a > - 1/2, 

Consequently the diffraction pattern should scale with DM ~ with 

- 1 > 6 > -1.5, rather than 5 = -1/2 as given in equation 6. Further 

measurements of a for high dispersion measure pulsars would be useful. 

Speeds of 50 krns” 1 are higher than values expected for the transverse 

components of earth and solar motions with respect to the local 

standard of rest, random hydrogen cloud velocities and differential 

galactic rotation over a few hundred parsecs. Some larger speeds 

o 

either in the medium (Wentzel 1970) or of the pulsar (Gunn and 
Ostriker 1970, Manchester e_t _al. 19 74) must contribute to the 
phenomenon . 
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Table 1 


Pulsar Observations Containing Interstellar Scintillation 


Experiment 

Goal 

Date Radio 

Frequency 

(MHz) 

Teles- Radiometer Sample 

cope Description Rate 

site- 
diam. (m) 

Time 

per 

Source 

(hr) 

Sources 

Collaborate 

Fluctuation 

1972 

Apr 

250-500 

GB-91 

250K 

transistor 

lp i 

0.7 

3 

None 

Flux density cI 

1972 

May \ 

250-500 

GB-43 

250K 

transistor 

60S ) 

5-8 

12 

J.R. Fishei 

Polarization j 



l 

1490-8085 

GB-23 

.12 OK 

paramps lOOOP^ 


1 

| J.M .Rankin 

Fluctuation 0 ' 

1971- 

-1974 

430 

AO-305 

1 50K 

pa ramp 

lp x 

0.5-3 

25 j 

| D . B . Campbe] 

Scattering 

1972 

Oct 

408 | 

GB-43 

| 250K 

transistor 

'1-5P 

5-12 

13 1 

JA.G. Lyne 






JB-76 





< 

|G.A. Zeissj 





1 

AO-305 







Scattering 

1973 

Feb 

111 

GB-91 

300K 

transistor 

lp i 

0.7 

5 | 

G . A . Zeissif 






AO-305 






N . R. Vandent 












T. A. Clark 

d 

Scattering 

1973 

| 

Jun 

834 1 

GB-43 

| 200K ; 

| paramps 

64P.J 

0.5 

1 

None 




Apr 

1667 


7 OK | 

i i 

i i 

128P 1 






i 

Jul 

4900) 


IIOOK 

! 1 

512P 1 ) 




Flux density 

1973 

May 

{ 250-500 

GB-43 

250K 

transistor 

60S ) 

5-8 

15 

J.R. Fislie: 





]l490-8085 

GB-26 

12 OK 

paramps lOOOP^j 




Scattering 

1974 

Feb 

1410 | 

( GB-43 

5 OK 

paramp ^ 

lp i 

0.3 

1 

M. Ewing 






mv-40 

100K 

par amp ( 






{ 



Table 1 (cont . ) 


Experiment 

Goal 

Date 

Radio Teles- 

Frequency cope 
(MHz) site - 

diam. (m) 

Radiometer 

Description 

Sample 

Rate 

Time 

per 

Source 

(hr) 

Sources 

Collaborators 

Position 

1974 May 

2695 

GB-26 

120K 

paramps 

0„ 3 

hr 

4 

5 

R. A. 

Sramek 

Flux 

density 

1974 May 

(775/968 
} 1440 

j GB-91 j 

i ) 

1 50K | 
50K \ 

paramps 

1 

4h 

hr 

0.1 

53 

J. R. 

Fisher 

a 

see 

Backer 

and Fisher 

(1974) 










b 

see 

Rankin , 

Campbell and Backer 

(1974) 









c 

see 

Backer . 

Rankin and 

Campbel 1 

(1974) 










d 


see Backer (1974) 


Table 2 

400 MHz Scintillation Spectra 


PSR 

P 

(s) 

DM 

(pc-cm~3) 

BW 

(MHz) 

NxM 

(P 1 ) 

m 

f e 

(mHz) 

t 

( loHz) 

f T 
e 

Data 

0031-07 

0.943 

10.9 

3.0 

1,024x1 

0.80 

1.3 

1 . 5 

1 

GBFF 

0301+19 

1.388 

15.7 

0.5 

8,192x1 

0.90 

0.88 

1.64 

10 

JBIS 

0329+54 

0.715 

26.8 

0.5 

16,384x3 

0.55 

1.48 

2.80 

19 

JBIS 

0525+21 

3.745 

50.8 

2.0 

1 ,024x1 

0.92 

12.5 

13.0 

52 

AOPN 

0540+23 

0.246 

72. 

0.1 

4 , 096x1 

0.35 

27.8 

52.0 

28 

AOPN 

0611+22 

0.335 

96.7 

0.1 

4.096x1 

0.30 

5.1 

7.0 

24 

AOPN 

0809+74 

1.292 

5.7 

3.0 

2,048x1 

1.00 

0.8 

3.0 

1 

GBFF 

0823+26 

0,531 

19 . 4 

0 . 5 

32,768x1 

0.95 

1.85 

4.70 

30 

GBIS 

0834+06 

1.274 

12.9 

2.0 

16,384x1 

0.75 

0,86 

2.00 

6 

JBIS 

0943+10 

1 . 098 

15.4 

2.0 

1,024x1 

0.55 

3.55 

5.30* 

4 

AOPN 

0950+08 

0.253 

3.0 

2.0 

8,192x5 

1.15 

0.19 

0.45 

2 

JBIS 

1133+16 

1.188 

4.8 

2.0 

8,192x1 

1 .07 

1.30 

2.60 

17 

JBIS 

1237+25 

1.382 

9.3 

2.0 

2,048x1 

1.00 

0.70 

1.75 

2 

AOPN 

1508+55 

0.740 

19.6 

0.5 

16,384x1 

0.30 

4.30 

5 . 00 3 

52 

JBIS 

1541+09 

0.748 

35.0 

2.0 

2,048x1 

£0.05 

<0.65 

- 

- 

AOPN 

1604-00 

0.422 

10.7 

2.0 

2,048x1 

o 

rH 

O 

VI 

*1.15 

- 

1 

AOPN 

1642-03 

0.388 

35 . 7 

0.5 

32,768x1 

0.60 

3.75 

7.00 

48 

JBIS 

1915+13 

0.195 

94. 

0.1 

2,048x1 

0.20 

5.00 

6.30 a 

2 

AOPN 

1919+21 

1.337 

12.4 

2.0 

8,192x1 

1.23 

0.64 

2.00 

7 

GBIS 

1929+10 

0.227 

3.2 

2.0 

16,384x1 

1.65 

1.07 

1.25 

4 

JBIS 

1933+16 

0.359 

158,5 

0.1 

1 ,024x1 

0.08 

8 . 20 

8.20 

3 

AOPN 

1944+17 

0.441 

16.3 

2.0 

2,048x1 

- 

*1.0 

- 

- 

AOPN 

2016+28 

0 . 558 

14.2 

0.5 

16,384x3 

0.75 

0. 58 

0.55 

16 

JBIS 

2020+28 

0.343 

24.6 

0.5 

2,048x1 

0.90 

2.85 

6.80 

3 

AOPN 


Table 2 (eont) 


PSR 

P 

(s) 

DM 

(pc-cm”3 ) 

BW 

(MHz) 

NxM 

( V 

m 

f e 

(mHz) 

ft 

( mHz ) 

f T 
e 

Data 

2021+51 

0.529 

22.6 

0.5 

16,384x1 

0.75 

0.69 

1.85 

6 

JBIS 

2045-16 

1.962 

11.5 

2.0 

4,096x1 

1.65 

(1.00) 

1 . 57 

11 

JBIS 

2217+47 

0.538 

43.5 

0.5 

16,384x3 

0.35 

2.12 

4 . 10 a 

56 

JBIS 

2303+30 

1 . 576 

49.9 

2.0 

2,048x1 

-^0.37 

(3.10) 

(3.0) 

1 

AOPN 

a 

Values 

may be 

limited by receiver 

noise rather 

than by 

intrinsic 

fluctuations 

as 


discussed in text. 




Table 3 



SCINTILLATION 

PARAMETERS AT 

2695 MHz 


PSR 

IX 

£ e 

(mHz) 

f T 
e 

Data 

0329+54 

0,38 

0.23 

4.5 

A 


0.90 

- 

- 

B 

0355+54 

0.38 a 

0.58 

4.5 

A 

0628-28 

0.69 

0.12 

2.5 

A 

0823+26 

0.94 

0.32 

8.5 

A 

0833-45 

0 . 10 a 

- 

- 

A 

0950+08 

^0,10 

- 

- 

A 


0.22 

- 

- 

B 

1133+16 

0.24 

0.24 

6 

A 


0.52 

- 

- 

B 

1237+25 

0.35 

0.36 

8 

A 

1706-16 

0.45 

0.35 

5.5 

A 

1749-28 

1 . 05 

0.14 

2 

A 

1929+10 

0.53 

0.18 

4 

A 


0.52 

- 

- 

B 

1933+16 

0 . 20 a 

0.25 

6 

A 

2016+28 

-0 . 67 

0.05 

<1 

A 

2021+51 

0.70 

0.10 

2.5 

A 


0.84 



B 


^'Suppressed from ~1.0 by 60~MHz bandwidth. 
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Fig. 1 
Fig. 2 

Fig. 3 


Fig . 
Fig. 

Fig. 


CAPTIONS 

Schematic spectrum of pulsar fluctuations with four 
regions (a-d) identified which are discussed in the 

text. 

Normalized variance-density spectrum for PSR 0329+54. 
Upper abscissa scale is for the high-resolution 
spectrum (solid) and lower abscissa scale is for the 
low-resolution spectrum (dashed) . A gaussian curve 
(dash-dotted) is drawn through the low-frequency 
end of the high-resolution spectrum. 

Normalized variance-denstiy spectrum for PSR 0823+26. 
Upper abscissa scale is for the high-resolution 
spectrum (solid) from 408-MHz data and lower abscissa 
scale is for the low-resolution spectrum (dashed) 
of 430-MHz data. A gaussian curve (dash-dotted) is 
drawn through the low-frequency end of the high- 
resolution spectrum . 

1 Normalized variance-density spectrum for PSR 1133+16. 
See caption to fig. 2 for further details. 

5 Normalized variance-density spectra for PSR’s 0329+54 
(a) , 0823+26 (b) , and 1133+16 (c) from 2695- MHz 
data. The scales are identical to those of figuies 

2 to 4, respectively. 

6 Scintillation bandwidths of pulsars near 400 MHz.. 
Upper limits are shown for three objects with 
arrowed points. The point outside the dotted 
lines at DM = 70 pc cm’ d is for the Vela pulsar. 



Fig. 7 


Fig. 8 


Scintillation indices of pulsars at 2695 MHz from 
observations of the author and from published 
data of Downs and Reichley (1971) 

Radio- frequency dependence of scintillation 
index estimates for nine pulsars. 
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APPENDIX I* 

VARIANCE-DENSITY SPifiCThA FOR 28 PULSARS 

In Hill and IV of this paper we discuss the shapes 
of the scintillation portion of normalized, variance- 
density spectra (VDS's) observed for pulsars near 400 MHz 
and relate the width of these functions to the dispersion 
measures, or column density of interstellar thermal electrons. 
In this prepublication document we present the VDS ’ s for 28 
pulsars upon which the discussion in § III and the measure- 
ments in §IV are based. 

For most objects a high-resolution VDFCsolid line and 
upper frequency scale) and a low-resolution VDS (dashed line 
and lower frequency scale) are given on a semi logarithmic 
grid in the following pages. The low-resolution frequencies 
are in units of cycles per pulse period (c/P ) . The pulsar 
designation and observing frequency (-ies) are noted in the 
upper righthand corner of each plot. When two frequencies are 
given, the first (408 MHz) refers to the high-resolution VDS 
and the second (430 MHz) refers to the low-resolution VDS . 

This division occurs since the 408-MHz data are of longer 
extent, allowing higher resolution and a more nearly station- 
ary VDS, while the 430-MHz data have superior sensitivity, 
^The - fol lowing pages will not be submitted for publication. 



revealing the short-term fluctuations in cases of weak signals. 
In these cases the high-resolution spectrum typically decays 
to a level of the system noise in the 408-MHz observations 
(e.g. , see PSR 0823+26). 

Error bars are given for both VDS ' s based on the number 
of degrees of freedom in the plotted points: two. for the 

sine and cosine components in the raw VDS . times _n for the 
smoothing performed before plotting. The error was taken 
to be a factor of (2 n) 1 ' 2 . For a stationary random signal, 
such as the typical system fluctuations, this error should 
give the s t a nda r d dev ia t ion o f t he points of a VDS f rom the 
(local) mean. In many cases we have not obtained sufficient 
data to provide a stationary estimate of the scintillation 
portions of the VDS ' s . This is discussed in more detail in 
§ § 1 1 1 and VIII. On the other hand, it is found that the low- 
resolution VDS’s do provide stationary estimates of the 
fluctuations intrinsic to the pulsar emission . Deviations 
in the low-resolution VDS's from a "white" continuum which 
are significant relative to the errors are related to 
periodic subpulse variations (e.g., see PSR 0823+26), or to 
bursts of emission (e.g. . see PSR 1944+17) ; these 



deviations are discussed in more detail by Backer, Rankin 
and Campbell (1974) . 

A third curve (dash-dotted line) is given which is 
a subjective fit of a half -gauss ian to the scintillation 
portion of the high- resolution VDS , when it can be identi- 
fied, as discussed in §111. In most cases the VDS follows 
a gaussian curve within the quoted statistical errors. A 
further constraint on the fitting was that the area of the 
gaussian be near unity. This fact explains why a poor 
"fit", is given for PSR 1919+21. It was felt that since 
unity for scintillation indices is appropriate in most 
objects, it is probably appropriate for all objects. The 
’’excess'’ variance-density at scintillation-like frequencies 
then would be due to variations intrinsic to t he pulsar . 

This hypothesis could be tested by simultaneously inspecting 
the correlation bandwidth of the fluctuations to determine 
what fraction of the fluctuations is due to the (narrow- 
bandwidth) scintillation. This experiment would be useful 
for all objects listed in table 2 with large ‘ est imated 


scintillation indices. 
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